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Abstract 

We investigate the "generalized Heron polynomial" that relates 
the squared area of an n-gon inscribed in a circle to the squares of its 
side lengths. For a (2m + l)-gon or (2m + 2)-gon, we express it as 
the defining polynomial of a certain variety derived from the variety 
of binary (2m — l)-forms having m — 1 double roots. Thus we ob- 
tain explicit formulas for the areas of cyclic heptagons and octagons, 
and illuminate some mysterious features of Robbins' formulas for the 
areas of cyclic pentagons and hexagons. We also introduce a com- 
panion family of polynomials that relate the squared area of an n-gon 
inscribed in a circle, one of whose sides is a diameter, to the squared 
lengths of the other sides. By similar algebraic techniques we obtain 
explicit formulas for these polynomials for all n < 7. 



1 Introduction 

Heron of Alexandria (c. 60 BC) is credited with the formula that relates the 
area A" of a triangle to its side lengths a, b, and c: 

K = a/ s(s — a)(s — b)(s — c) 

where s = (a + b + c)/2 is the semiperimeter. For polygons with more than 
three sides, the side lengths do not in general determine the area, but they do 
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if the polygon is convex and cyclic (inscribed in a circle). Brahmagupta, in 
the seventh century, gave the analogous formula for a convex cyclic quadri- 
lateral with side lengths a, b, c, and d: 

K = y/{s-a){s-b){s-c){s-d) 

where s = (a + b + c + d)/2. See [2] for an elementary proof. 

Robbins [?| found a way to generalize these formulas. First, drop the 
requirement of convexity and consider the square of the (signed) area K of 
a possibly self-intersecting oriented cyclic polygon. For this purpose we can 
define the area enclosed by a closed curve C to be § c xdy. Second, express 
the relation between K 2 and the side lengths as a polynomial equation with 
integer coefficients. Given a cyclic polygon, one can permute its edges within 
its circumscribed circle without changing its area, so the polynomial will be 
symmetric in the side lengths, and in fact it can be written in terms of 1QK 2 
and the elementary symmetric functions <7j in the squares of the side lengths. 
For instance, the Heron and Brahmagupta formulas can be written 

16K 2 - 4o- 2 + of ~ e ■ 8^ = 

in which e is for a triangle, 1 for a convex quadrilateral, and —1 for a non- 
convex quadrilateral. Hence all cyclic quadrilaterals satisfy the polynomial 
equation (16K 2 — 4er 2 + af) 2 — 64cr 4 = 0. The general result is the following. 

Theorem 1. [7] For each n > 3, there is a unique (up to sign) irreducible 
polynomial a n with integer coefficients, homogeneous in n + 1 variables with 
the first variable having degree 2 and the rest having degree 1, such that 
a n (16K 2 , a\ , . . . , a%) — whenever a\, . . . , a n are the side lengths of a cyclic 
n-gon and K is its area. 

The polynomials a n are now known in the literature as generalized Heron 
polynomials. For certain sets of n side lengths, as shown in jjj, one can find 
up to A n distinct squared areas, where 

n fn — 1\ n _ 2 

Hence one expects that a n has degree A n in its first variable. This conjecture 
of Robbins, and two others made in [2j , have recently been established. We 
summarize them in Theorem 121 
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Theorem 2. The polynomial a n is monic in 16K 2 and has total degree 2A n . 
If n is even, then a n = /3 n /3* where f3 n is a polynomial in 16K 2 , a%, . . . , o n -\, 
y/dn, where y/dn = ai---a„ and /3* is (3 n with ^fo^ replaced by —y/dn- 

See (HI or §5] for the degree, [S] or [T] for monicity, and [S] for the factor- 
ization when n is even. Robbins' main interest, however, and the motivation 
for our research, was to find reasonably explicit formulas for all a n and f3 n . 

In [7j, Robbins found formulas for and (3$ that have a curious form. To 
present them concisely, we introduce the crossing parity e of a cyclic n-gon. 
Assume that the n-gon has vertices v±, . . . , v n in the complex plane and 
circumcenter 0. For odd n let e = 0, and for even n let e = ±1 be 

e = sign[(l - f 2 /fi)(l - vz/v 2 ) • • • (1 - vi/v n )], 

which is well-defined because the product is real and nonzero if all edges 
have positive length. (This definition agrees with the previous definition of 
e for n e {3, 4}.) Now assume n G {5, 6}. Define u 2 = —4K 2 , and make the 
substitutions 

ti = ox, 

t 2 = -o- 2 + |*i - u 2 , 

t 3 = Or 3 + \t X t 2 - t ■ 2^, (1) 
1 2 I 

ti = —0-4 + jt 2 + e • Ti-y/cre, 

Then, for any cyclic pentagon or hexagon of the given crossing parity, the 
cubic polynomial u 2 + t%z + t^z 2 + t^,z 3 has a double root, so its discriminant 
vanishes: 

t\t\ - Au 2 t 3 4 - At% + I8u 2 t 3 t 4 t 5 - Wu\t\ = 0. 

When the are expanded, this discriminant is a polynomial of degree A 5 = 7 
in u 2 , and hence in 16K 2 . Multiplying it by 2 18 makes it monic in 16K 2 and 
yields 05, (3q, or /5| according to whether e is 0, +1, or —1. 

In JJJ we generalize this construction. Fix n and the crossing parity e, 
and let m = \_(n — 1)/2J . We introduce auxiliary quantities u 2 , u m , 
with u 2 = —4K 2 , and inductively define certain polynomial expressions £j in 
the Oj and Uj with j < i. For n = 5 or 6, these definitions reduce to (JTJ). 
Corollary El then says that the polynomial 

Pn{z) — U 2 + • • • + U m Z + t m+ iZ + • • • + t 2m+ \Z 
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is divisible by the square of a polynomial of degree m — 1. In other words, for 
any values of the and Uj coming from a cyclic n-gon, P n {z) has m — 1 double 
roots over C (counting with multiplicity, and including roots at infinity). 
Such polynomials form a variety of codimension m — 1, defined locally by 
m—1 equations. So, if we regard 113 through indeterminates and 

expand each ti in terms of the <7j and Uj, we can in principle eliminate the 
m — 1 unwanted quantities %, . . . , u m and recover a single polynomial, which 
is a2m+i, /?2m+2, or /?2m+2 depending on e. In HJwe carry out this program for 
m = 3 to obtain formulas for a-j and as, the generalized Heron polynomials 
for cyclic heptagons and octagons. 

There is another family of area polynomials susceptible to the same anal- 
ysis. Call an (n+ l)-gon semicyclic if it is inscribed in a circle with one of its 
sides being a diameter. Its squared area satisfies a polynomial relation with 
the squares of the lengths of the other n sides; the degree in the squared area 
turns out to be 



Theorem 3. For each n > 2, there exists a unique monic irreducible poly- 
nomial a' n with integer coefficients, homogeneous in n + 1 variables with 
the hrst variable having degree 2 and the rest having degree 1, such that 
) = whenever the lengths of the sides of 

a semicyclic (n + l)-gon excluding a diameter, and K is its area. The total 
degree of a' n is 2A' n . 

The proof that a' n exists and is unique (without assuming monicity) fol- 
lows the proof of Theorem in |7j almost verbatim, and the argument in £Q 
shows that a' n is monic. We establish the degree by an elementary argument 
in £0 which is independent of the rest of this paper. 

Cyclic and semicyclic polygons are similar in many ways. For instance, 
just as the polygon of largest area one can make with n given side lengths 
is convex and cyclic, the polygon of largest area one can make with n given 
side lengths and one free side is convex and semicyclic. We will adduce many 
algebraic similarities in the following sections. For now we just observe that 
the polynomial ctg, which can be worked out by hand, also takes the form of 
a discriminant: if U2 = —AK 2 , then 




«3 = 16 discr 2 (z 3 + a x z 2 + (a 2 + u 2 )z + cr 3 ) . 



(2) 
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2 The Main Identity 



All our area formulas are based on a generating function identity that relates 
the symmetric functions (7j in the squared side lengths to certain quantities 
Tj that arise in Robbins' proofs of the pentagon and hexagon formulas. The 
identity, Theorem holds for both cyclic and semicyclic polygons and for 
both odd and even n. 

Suppose we have a cyclic n-gon or semicyclic (n + l)-gon inscribed in a 
circle of radius r centered at the origin in the complex plane. Let its vertices 
be Vi, . . . , v n and v n+ i = Svi, where 6 = 1 for a cyclic n-gon and 5 = — 1 
for a semicyclic (n + l)-gon. Introduce the vertex quotients q^ = v i+ i/vi for 
i — 1, . . . , n, and let r , r 1; . . . , r n be the elementary symmetric functions 
of the qi. Then r = 1 and r n = qiq 2 ■ ■ - q n = 5. Elementary geometry yields 
the equations 

a 2 = r 2 (2-q l -qr 1 ), 1 < i < n, (3) 

16K 2 = -r\ gi + ... + q n -q-' q^f (4) 

= -r 4 (r! - 5r n ^f. 

Using © one can express each <jj in terms of r and the Tj. Let g(y) = 
y 2 + (x/r 2 — 2)y + 1. Observe that x is one of the values a 2 exactly when g(y) 
has one of the vertex quotients root, or in other words, when g(y) has 

a common root with the polynomial f(y) = YYi=i(y ~ 1i) = Sr=o( — ^) %T iV n ~ l - 
Hence the resultant of f(y) and g(y) is a constant times 

n n 

h{x) = ]j(x - a 2 ) = Y^{~ir^x n -\ 

i=l i=0 

and the coefficient of x n reveals that the constant is 5r~ 2n . By expanding 
the resultant, one finds that each Oi is r 2% times a quadratic polynomial in 
the Tj. A particularly simple example is 

a n = 5(-l) n r 2n (r - n + r 2 - - • • ± r n ) 2 . 

If n is even and 5 = 1, then y/dn is expressible in terms of r 2 , the Tj, and the 
crossing parity e: 

y/on = \vi - v 2 \ ■ ■ ■ \v n - v n+1 \ = r n |1 - gi| • • • |1 - q n \ 

= r n e(l -q x ) ■■■{!- q n ) = r"e(r - n + r 2 - • • • + r n ). ( ''' 
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So far we are following |7j except for the addition of the semicyclic case. 

Consider now the involution that reflects the polygon in the real axis. 
This operation preserves the squared area and the side lengths, but it re- 
places each qi with q~ l and hence replaces each r, with <5r n „j. Because each 
<Tj is a quadratic form in the Tj preserved by the involution, it can be uniquely 
decomposed into a two parts: a quadratic form in symmetric linear combi- 
nations of the Tj, and a quadratic form in antisymmetric linear combinations 
of the Tj. When we perform this decomposition on the whole generating 
function x Y a ii each part factors in a surprising way, which our main 

identity records. 

To write the identity explicitly, we need the following linear combinations 
of the Tj, for < k < n/2: 



* = £(-!)'( 

k 

e* = £(-!)* 

i=0 



n - 1 - 2k + i 



n-2k + i\ (n - 2k + % - 1 

i M i-i 



(6) 

ijk-i + T n - k+i ). (7) 



Let D(x) = diX % and E(x) = e i x 



Theorem 4 (Main Identity). For a cyclic n-gon or semicyclic (n + l)-gon 
of radius r, with 5 = 1 or — 1 respectively, the symmetric functions Oi of the 
squared side lengths and Ti of the vertex quotients are related by 



-X (Ti 



\E{r 2 xf + (r 2 x - \)D{r 2 xf. 



i=0 



Proof. When the o~i are expanded in terms of the Tj, both sides of the main 
identity become polynomials in r 2 x, so we may assume r = 1. The left-hand 
side is then 

n 

5 ■ ^(-x)V, = 5x n h{x- 1 ) = x n Res(/, g) 

i=0 

where f(y) = Er=o( -1 ) V ^ n ~ t and d(y) = (y ~ I) 2 + 

We calculate the resultant using its PGL(2)-invariance and other stan- 
dard properties j3]. Make the change of variable y — (z — l)/(z + 1) so that 
the roots of g are related by z i— > —z instead of y i— > y^ 1 . We obtain the 
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polynomials 

n 

f*(z) = (z + l) n f (*=i) = ^(-l)V,(z - + 1)*, 

i=0 

g*(z) = (z + l)*g(£)=x- 1 {* l + 4x-l), 
and the transformation has determinant 2, so 

x n Res(/, = 2- 2ri x" Res(/*, #*) = 2~ 2n /* ( V / T^4i) /* ( -vf^Sc) . 
Write f*(z) = fo(z 2 ) + zfi(z 2 ), separating even and odd powers of z. Then 

f (l-Ax) 2 -(l-Ax)f 1 (l-Axf , 

which explains the form of the main identity. It remains to evaluate / (1 — 4x) 
and — Ax). We consider only / , as f\ is similar but simpler. 

It helps to introduce the Fibonacci polynomials F n {x) = Ylt^o^ ( n 7*) x *' 
which count compositions of n by l's and 2's. They satisfy the recurrence 

F n (x) = F n _i(x) + xF n _ 2 (x), n > 1, (8) 

with Fq(x) = 1 and F n (x) = for n < 0, and have generating function 

F{x; t) = J2 F n{x)t n = (1 - t - xt 2 )- 1 . 

n 

The generating functions for the c4 and can then be written 

Ln/2j 

D ( x ) = E ( T i ~ T n -i)x l F n _2i-i{-x), 
i=0 
\n/2\ 

E ( x ) = E (n + T n _i)x l (F n _ 2i (-x) - xF n _ 2i _ 2 (-x)) . 

i=0 

To evaluate / (1 — Ax), first rewrite fo(z 2 ) = \ {f*{z) + f*(—z)) in terms 
of the sums T; L + r n _j. If we let 6 = | if 2i = n and = 1 otherwise, then 

[n/2j 

/o(^ 2 ) = E (-i) B-i *fo + - E 

i=0 i 

7 



x n Res(/,( 7 ) = 2- 



(n-2i\ . 



Next, we want to substitute z 2 = l — 4x and evaluate the sum over j. Writing 
m = n — 2i, we can compute the generating function 

E r E (i - *Y - = f 1 - 2 <> 

m>0 j v J/ 

by interchanging sums and simplifying. Thus we get 

Ln/2j 

/o(l-4x) = ^ (-l) i+n 9(r t + r n ^)(4xY2 n - 21 (F n _ M (-x) - ^^(-x)) . 

8=0 

The recurrence (jHJ shows that 6 1 (F„_ 2 j(— x) — ~_F n _2i_i(— a?)) is equivalent 
to | (F n - 2i (-x) - xF n -2i-2(—x)), and so 

f (l-4x) = (-l) n 2 n ~ 1 E(x). 

Likewise /i(l — 4x) = (— l) n 2 n_1 D(:c), and the identity follows. □ 



3 Consequences of the Main Identity 

The main identity tells us how to generalize the definition of the quantities ti 
and Uj that were so useful in simplifying the pentagon and hexagon formulas. 
Cyclic n-gons have d = r — r n = and e = t + r n = 2 by © and (JZj) , so 
the expansions of E(r 2 x) 2 and D(r 2 x) 2 include linear terms in the but not 
in the dk- The substitutions that replace the <7j with the £, will first isolate 
and then eliminate the variables e^, and leave us with expressions relating 
the t{ and Uj to the radius r and the variables dk- The main identity will then 
express the algebraic relationship among the ti and Uj as the factorization of 
a single polynomial P n (z). 

Corollary 5. Given a cyclic n-gon of crossing parity e and radius r, let 
m = [(n — 1)/2J and let Uj = r 2 - 7 ^i=i(l^i — di-i)dj-i f° r j > 1- Inductively 
define t = —2 and 

for j = 1, . . . , 2m + 1. Then tj = — e^r 2 - 7 for < j < m, and the polynomial 

P n (z) =u 2 + u 3 z + --- + u m z m ~ 2 + t m+1 z m - 1 + ■■■+ t 2m + x z 2m ~ x 
factors as (| — r 2 z)[z~ l D{r 2 z)] 2 . 
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By (JHJ) and (gj), we have u 2 = \r^d\ = |r 4 (ri - r n _i) 2 = -4K 2 , and 
u\ = by definition. Thus the tj and Uj in Corollary 03 agree with those 
defined in ^ 

Proof of Corollary^ We have e = t + r n = 2 by the definition (J7|), so 
t = — eor°. Now let 1 < j < m. The coefficient of x J in ^E(r 2 x) 2 is 

i=0 ^ i=l ^ 

by induction on j. The coefficient of x- 7 in (r 2 x — j)D(r 2 x) 2 is — u,, so the 
equation tj = follows by comparing coefficients of x J in the main 

identity 

For j > m we must consider the coefficient r 2m+2 e m+ i of x m+l in E{r 2 x). 
If n = 2m + 1, then E has degree m by definition so this coefficient is zero. 
But if n = 2m + 2, then the coefficient is 



r 



(m+l \ 
2r m+1 + ^(-l) J 2(r m+1 _, + r m+1+i )j 



i=l 

\m+l 



= 2(-l) m ^e^ 

by (0) and (0). So for m < j < 2m + 1, the coefficient of x J in jE(r 2 x) 2 is 

m m 



r 

i=j—m i=j—m 

and this holds whether n is odd or even because e = when n is odd. Thus, 
by the main identity, — tj is the coefficient of x^ in {r 2 x — ^)D{r 2 x) 2 for 
j = m + 1, . . . , 2m + 1. 

We now see that {r 2 x — ^)D{r 2 x) 2 , a polynomial of degree 2m + 1 whose 
two lowest terms vanish, is exactly — x 2 P n (x). □ 

There is a geometric argument that Corollary El contains enough infor- 
mation to recover a n . To simplify the explanation, assume n = 2m + 1 and 
m > 2. The nonzero polynomials of degree up to 2m — 1 that have a squared 
factor of degree m — 1 naturally form a projective variety of codimension 
m — 1 in p 2m - 1 ; which is irreducible because it is the image of P 1 x p m_1 
under a regular map. Hence the affine variety X m C A 2m of such polynomials 
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(now including the zero polynomial), which has the same ideal, is also irre- 
ducible. The substitutions that write t m+ \, . . . , t 2m +i in terms of the <7; and 
Uj amount to a morphism / : A 3m — > A 2m , which is a product bundle with 
fiber A m . This is because, for any point (m 2 , . . . , u m , t m+ i, . . . , tzm+i) in the 
range and any given values of a%, . . . , a m , the values of cr m+ i, . . . , cr 2m+ i are 
uniquely determined as polynomial functions of the other variables. Hence 
/ -1 (Jf m ) pd X m x A m is irreducible and of codimension m—1. Finally, when 
we apply the projection 7r : A 3m — > A 2m+2 that eliminates the m — 2 variables 
«3, . . . , u m , the closure of the image 7r(/ _1 (X m )) is an irreducible variety of 
codimension at least 1 that contains V(a n ), so it must equal V(a n ). 

Corollary therefore reduces the problem of finding a n to two subprob- 
lems: finding the defining equations of the variety X m , and then, after ex- 
panding t m +i, . . . , t2m+i in terms of the <j{ and Uj, eliminating the m — 2 
variables 113, . . . , u m . 

The application of the main identity to semicyclic polygons is similar but 
slightly different. In this case e = t + r n = and d = r — r n = 2, so the 
main identity involves linear terms in the c4 but not the e^. Our definitions 
of ti and Uj are therefore designed to extract and eliminate the variables d^. 
Again we can distill the relationship among the ti and Uj to the factorization 
of a polynomial -P^(z). This time, due to the factor {r 2 x — |) in the main 
identity, the expression for ti explicitly includes r 2 , so there remains one more 
unwanted variable to eliminate for a given n. 

Corollary 6. Given a semicyclic (n + l)-gon of radius r, let m = \_(n — l)/2j 
and let Uj = r 2 i YliZi e i e j-i/^ f° r 1 < j < m - Inductively define to = —2 
and 

fi = (-i)' + V i+ £ 

l<i<m 
l<j—i<m 

for j = 1, . . . ,n. Then tj = —djr 2 ^ for < j < m, and the polynomial 

P^z) =u 2 + u 3 z + --- + u m z m ~ 2 + t m+l z m ~ l + ■■■ + t n z n - 2 
is the square of E{r 2 z)/2z. In particular, t n = if n is odd. 

Once again u 2 = \r A e\ = \r A {ji + r n _i) 2 = — AK 2 by (jH), since now 
5 = -l. 



— r 



E 



-U 



J * 



0<«-l<m 
0<j~i<m 



if j < m, 
if j > m, 
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Proof. As in Corollary |3J the claims follow from equating coefficients of x J 
in the main identity and inducting on j to evaluate tj for < j < m. If 
n = 2m + 1, the degree of E(x) is just m, so t n = 0. □ 

The polynomial P^{z) contains m — 1 unwanted variables, namely r 2 and 
ii 3 , . . . , u m . If n = 2m + 1, then P^z) is a polynomial of degree 2m — 2 that 
is a square, which gives rise to m — 1 equations in its coefficients, and we 
have the additional equation t n = 0. If n = 2m + 2, then P^(^) is a square 
of degree 2m, which yields m equations. In either case Corollary |B] holds 
enough information, in principle, to derive the area formula a' n . As before, 
one can make this claim precise using some algebraic geometry. 

4 Explicit Formulas 

In this section we apply the results of §21 to produce area formulas for cyclic 
heptagons and octagons, and also semicyclic quadrilaterals, pentagons, hexa- 
gons, and heptagons. 

Because the degree of the generalized Heron polynomial a n is exponential 
in n, and the number of terms could be exponential in n 2 , there is some ques- 
tion as to what constitutes an explicit formula. Our formulas have concise 
descriptions, and if a polygon is given with exact (for instance, rational) side 
lengths, the polynomial satisfied by its area can be computed exactly using 
standard operations such as evaluating determinants. 

First let us apply Corollary [S] to the cases n G {7, 8}. It gives us a binary 
quintic form 

x 5 P n (y/x) = u 2 x 5 + UzX A y + t A x 3 y 2 + t 5 x 2 y 3 + t Q xy A + t 7 y 5 

whose coefficients are polynomials in u 2 , u 3 , <7i, . . . , 07 and perhaps 
and which, when its coefficients are evaluated for any cyclic ra-gon, has two 
linear factors over C of multiplicity two. The condition for a quintic form Q 
to factor in this way is given by the vanishing of a certain covariant C, which 
in the notation of transvectants jH] is 

C = 2Q(H,t)^ + 25H(Q,t)^ + 6 Qz 2 , H = (Q,Q) {2 \ t = (Q 1 Q)^. 

Here (f,g) {d) = ^ =Q (-l) i (f)(d d f/dx i dy d - i )(d d g/dx d ~ i dy i ). This fact about 
quintics is presumably classical, but we have not yet found a reference. In 
any case, C is a form of degree 9 in {x, y} whose coefficients are forms of 
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degree 5 in the coefficients of the original quintic, so its coefficients give us 
ten degree-5 polynomials in u 2 , Us, t±, t 5 , t 6 , t 7 that must vanish. These same 
ten polynomials can be obtained as the Grobner basis, with a graded term 
ordering, for the ideal of the variety of quintic forms that factor as a linear 
form times the square of a quadratic. 

To obtain the desired relation between u 2 and the cr.j, we must expand 
the coefficients of C as polynomials in w 3 and then eliminate %. We can do 
this most explicitly using resultants with respect to u%. The two simplest 
coefficients of C are 

F = u\t\ - Au 2 t\ - 4^3*4*5 + l8u 2 U 3 t\t b - 27^2*4*5 

+ (8W3 - A2u 2 ulU + "i<ou\t\ + hAulush - 80uij* 6 )* B 
+ {8u 2 u\ - 30ulu 3 t 4 + 50ult 5 )t 7 , 

of total degree 18, and 

G = u\t% - Au 2 t% - 4u 3 3 t 2 5 + 18u 2 u 3 t 4 tl - 21u\t\ 

+ (2u 3 t4 — 8u 2 u 3 t 4 — 6u 2 u 3 t 5 + 36u 2 t 4 t 5 — 8u 2 u 3 te)te 
+ (I6M3 - 7Au 2 ujU + tt$u\t\ + 110u|u 3 t 6 - 200u|t 6 )*7, 

of total degree 19. Let P *—>■ P denote the operation of expanding the £j in 
terms of 1*2, M3, and 01, . . . , as specified by Corollary El This operation 
preserves total degree. Both F and G have degree 6 in U3. Their resultant 
with respect to 113 therefore has total degree 6 x 19 = 114, and it must have 
the polynomial a 7 of total degree 2A 7 = 76 as a factor. 

The resultant Res(F, G) seems to be too large to compute and factor 
explicity, but we can describe its unwanted factors as follows with a little 
computer assistance. First observe that every term in F and G is divisible 
by either u 2 or U3, and hence the same is true of F and G. It follows that 
Res(F, G) is divisible by u 2 . In fact u\ | Res(F, G), as we will see in Lemma[7| 
below. Next, consider the polynomials 

Fx = 4«3 - 15u 2 u 3 t 4 + 25^5, 

Gx = 7u\H/x — 20u 2 tl — 5u 2 u 3 t 5 + lOOu^te, 

which are closely related to the coefficients of t 7 in F and G. Specifically, 
Fx =j2u 2 )- 1 dF/dt 7 and^Gi = u 2 1 (2u 3 F 1 - dG/dt 7 ). We will show that 
Res(Fi,Gi) divides Res(F, G). 
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First we claim that if F l = G l = 0, then F = G = 0. The ideal (F 1 , G t ) 
does not contain F and G, but by some easy calculations, it does contain 
U2F, U3F, U2G, and U3G. If F\ = G\ = 0, then all four of these polynomials 
vanish; so if either tt 2 ^ or M3 ^ 0, we must have F = G = 0, while 
if u 2 = u 3 = 0, we already know that F = G = 0^ This establishes the 
claim. It follows that F ± — G\ — implies F = G = 0. Consequently 
wherever Res(F 1 ,Gi) vanishes, so does Res(F, G). Algebraically, this means 
that every irreducible factor of Res(Fi, G±) divides Res(F, G). The resultant 
of F\ and G\ with respect to M3 is simple enough to compute explicitly. It 
has total degree 30, and it factors as u 2 times an irreducible polynomial in 
Q[u 2 , <7i, . . . , 07] of total degree 24. 

Thus, not only does Res(F!,Gi) divide Res(F, G), but u\ Res(Fx, G\) 
does also. The quotient by the latter polynomial has total degree 114 — 8 — 
30 = 76 = 2A7, so it must be a scalar multiple of the desired polynomial «7, 
Ps, or there are no more unwanted factors. The scalar can be computed 
by setting <j 2 , ■ ■ ■ , 07 to zero, and we find that 

2 101 5 5 Res(F,G> 3 ) 
u\ Res(Fi, Gi, M3) 

is 0:7, [3$, or /3g, according to whether the crossing parity e is 0, +1, or — 1. 
It remains only to prove the following lemma. 

Lemma 7. With the definitions above, u 7 2 \ Res(F, G, -u 3 ). 

Sketch of proof. By direct calculation on a computer, u 7 2 divides Res(F, G, 113) 
but u 2 does not. The only component of V(F, G) lying on the hyperplane 
u 2 = is the linear variety V(u 2 ,u 3 ), so V(F) and V{G) must intersect with 
multiplicity 7 along V(u 2 ,u 3 ) C A 6 . Now, assuming n = 7 for definiteness, 
pull back via the projection 71 : A 9 — > A 6 that maps (-u 2 , m 3 , oi, . . . , (j 7 ) 1— > 
(m 2 , m 3 , t 4 , . . . , ^7). Because 71 is smooth, the intersection multiplicity of V(F) 
and V(G) along 7r _1 \/(-u 2 , M3) = V r (-u 2 ,'U3) C A 9 is also 7. For fixed generic 
values of 01, . . . , 07, we therefore have u\ \ Res(F, G, u 3 ). We conclude that 
this divisibility holds globally as well. □ 

For the rest of this section, we turn our attention to semicyclic (n+l)-gons 
with n — 3, 4, 5, and 6. To state the area formulas most cleanly we introduce 
a notion of parity for semicyclic polygons. Let n be even, and observe that the 
quantities e 1 =T 1 + r n _i = E % - E and \ e n/2 = ZX -1 )* 7 "* = IK 1 _ ft) 
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are both pure imaginary. (Compute their complex conjugates using q i = q i 1 .) 
Hence their product is real. Let e G {—1, 0, +1} be its sign. Then we have 

e \K\ y/an = e-\r 2 |n + r n _i| • |«i - f 2 | |u 2 - «s| • • • |u„ - w i| 
= e • \r 2 |n + r n _i| • r n |1 - g x | ■ ■ • |1 - q n \ 
= r n+2 ■ i ei • \e n/2 . 

Define w = 2e\K\ y/dn = e\/u 2 t n for n even, and let w — for n odd. Our 
formulas for a' 4 and «g will factor when written in terms of w rather than 
<j n . We do not know whether this type of factorization occurs in general. 

For n G {3,4}, it is simplest to use the main identity directly. Defining 
C2 = if n = 3, we have 

\E{r 2 xf = ye\x 2 + y ei e 2 x 3 + ±r 8 e|r 4 
= u 2 x 2 + 2wx 3 + ■ ■ • 

so, by Theorem HI the cubic 1 — <J\X + (cr 2 + u 2 )x 2 — (cr 3 — 2w)x i factors 
as — (r 2 x — j)D{r 2 x) 2 . In particular, its discriminant vanishes. Replacing x 
by —x^ 1 , we recover equation (J2J) for n = 3, and for n = 4we have factored 
a' A as the product of two discriminants /3' 4 and ((3' 4 )* corresponding to e = +1 
and e = — 1 respectively. 

For larger n we need Corollary El For n — 5, it says that M2 + £3*2 + 
z 3 is the square of the linear polynomial E(r 2 z)/(2z), which yields 
the two equations t\ — 4^4 = and £5 = 0. Their degrees in r 2 are 6 and 5 
respectively, so their resultant with respect to r 2 has the correct total degree 
2A' 5 = 30. (Remember that r 2 has degree 1.) It remains only to scale the 
resultant to be monic in —4u 2 , and we get 

a' 5 = \ Res(tg - 4u 2 t 4 , t 5 , r 2 ). 

For n — 6, Corollary |U1 gives us the factorization 

u 2 + t 3 z + t 4 z 2 + t 5 z 3 + t 6 z 4 = |r 4 (ei + e 2 r 2 z + e 3 rV) 2 . 

Using w = ir 8 eie 3 , we derive the equations 

u 2 t 5 - t 3 w = 0, 
u 2 + + (^4 - 2w)z 2 = |r 4 (ei + e 2 r 2 z) 2 , 
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the second of which implies that t 2 — 4u 2 (t4 — 2w) = 0. Thus we can form 
the resultant of t 2 — 4-u 2 (^4 — 2w) and u 2 t 5 — t 3 w to eliminate r 2 and obtain 
a multiple of the desired area formula. The resultant is small enough to 
compute and factor symbolically, and we obtain a' 6 = (/3q) * where 

, _ Resjtl - Au 2 {U - 2^u 2 t 6 ), u 2 t 5 - t 3y /u 2 t 6 , r 2 ) 
^ ~ 4u° 

and (/3g)* is f3 r 6 with the opposite sign on y/u 2 t e . 



5 Degree Calculations 

In this section we show by elementary means that the homogeneous polyno- 
mials a n and a' n have degrees 2A n and 2A' n respectively. First we explain why 
the degrees cannot be smaller. In [7j, Robbins shows that deg(a n ) > 2A n by 
constructing A n cyclic n-gons with generically different squared areas from 
a given set of edge lengths. He takes the edge lengths to be nearly equal if n 
is odd, and adds a much shorter edge if n is even. For semicyclic polygons, 
we can take the edge lengths to be nearly equal if n is even; the argument of 
[7j then yields the desired number of semicyclic n-gons. 

Suppose now that n is odd. It is not necessary (and in fact not possible) to 
construct inequivalent semicyclic polygons with given positive real edge 
lengths dj. It suffices instead to construct A n configurations (r, qi, . . . , q n ) of 
complex numbers satisfying 

a 2 j =r 2 {2-q j -qr 1 ), j = l,...,n; (9) 

and qi ■ ■ ■ q n = —1, since it is from these equations, together with the relation 
16K 2 = —r^i^qj — ^qj 1 ) 2 , that one derives the existence and uniqueness 
of the irreducible polynomial a' n . In our configurations r is always real and 
positive, but sometimes r < minlaj/2}, in which case the qj are negative 
real numbers instead of complex numbers of norm 1. The plan is to regard 
each qj as a function of r by choosing a branch of equation @, and then find 
values of r such that q\ ■ ■ ■ q n = —1. 

Let n = 2m + 1, let the first 2m edge lengths be large and nearly equal, 
and let a n = 2. To find configurations with r > max{aj/2}, choose arbitrarily 
whether < argg n < n (the short edge goes "forward") or — n < argg„ < 
("backward"), and likewise choose a set of k < m of the long edges to 
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go backward. Then there exist m — k semicyclic polygons with the given 
edge lengths and edge directions whose angle sums ar g<?j are 71 , 37r, . . . , 
(2m — 2k — 1)tt. (Apply the Intermediate Value Theorem to J^arggj as r 
varies from maxja.,/2} to oo.) The total number of such configurations is 

fc=0 v 7 V 7 

To find configurations with r < min{aj/2} = 1, choose the branch qj < 
— 1 for exactly m of the long edges, and choose the branch qj > —1 for the 
other long edges. Let £j — +1 or e,- = — 1 respectively. As r — > 0, the product 
Ql ■ • ■ Q2m approaches the constant Y\j=i ^ > an d hence q\ - ■ • q n approaches 
if q n > —1 or — oo if q n < —1. By choosing the branch for q n according to 
whether q± - ■ ■ q2m, evaluated at r = 1, is greater or less than 1, we guarantee 
that qi ■ ■ ■ q n = — 1 for some intermediate value of r. Thus we obtain another 
|( 2m ) configurations. (The factor of 1/2 is present because inverting every 
qj preserves the radius and the squared area; it corresponds to reversing the 
orientation.) The total number of configurations is therefore at least 

. , . [2m\ n ( n — l\ . , 

To establish matching upper bounds on the degrees of a n and a' n , we 
proceed indirectly. First we revive an argument of Mobius from the 19th 
century which produces a polynomial of degree A n that relates r 2 for a 
cyclic polygon to the squared side lengths. Hence there are generically at 
most A n circumradii for a given set of edge lengths. For generic side lengths 
(in particular, no two equal) and a radius r that admits a solution (gi, . . . , q n ) 
to the system of equations (JHj) and q\ • • • q n = 1, the solution is unique up to 
inverting all the qj. (Any other solution would differ by inverting a proper 
subset of the qj, so those qj would need to have product ±1.) Thus, because 
r and the qj determine the area, there are generically at most 2A n possible 
signed areas, so deg(a n ) < 2A n . The same argument applied to semicyclic 
polygons will yield deg(a^) < 2A^. 

Given a cyclic polygon with circumradius r and side lengths 2y 1 , . . . , 2y n , 
let 6j = sin^ 1 (yj/r) be half the angle subtended by the j th side. Then, for 
some choice of signs 62, ■ ■ ■ e n (namely, 6j is +1 or —1 according to whether 
the j th side goes "forward" or "backward" relative to the first side), the sum 
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6i + e 2 02 + • • • + e n 6 n is a multiple of it. Therefore 

n r " sin (^ + +■■■+ tj n ) = o. (io) 

The factors of r make this a polynomial relation over Q between r 2 and the 
squared side lengths. To see why, introduce the variables Xj = r cos 9j = 
(r 2 - yf) 1/2 and rewrite (JTUJ) as 



Hi 



j=i i=i 



(11) 



where E\ — 1. The left-hand side of (jllj) has a great deal of symmetry. 
Obviously, flipping the sign of z/j is equivalent to negating €j. Flipping the 
sign of any Xj is equivalent to flipping 6j and negating each product over j. If 
j = 1, we can restore the condition e\ = 1 by flipping every 6j and negating 
every bracket. All these operations just permute and possibly negate all 
the 2 n_1 bracketed factors, so they leave the overall expression unchanged. 
Therefore, in the expansion of (fTTjl . each Xj occurs only to even powers, and 
hence each x 2 can be replaced by r 2 — y 2 . Likewise each yj occurs only to 
even powers. Thus we obtain a polynomial equation M(r 2 , yf,..., y 2 ) = 0. 

The remaining part of Mobius' argument uses series expansion to find the 
degrees of the leading and trailing terms of M. Fix the €j, and rewrite the 
bracketed factor of (fTTj) as 

n n 

II {\l r2 -yl + ie iVi) - II {\l r2 -y 2 3 - ie oVi 

j=l 3=1 

To find the term of highest degree in r, expand around r = oo; the highest 
terms cancel, so the degree is n — 1. To find the term of lowest degree, expand 
around r = to get 

" / r 2 \ n ( r 2 \ 

Its initial term has degree mm(k,n — k) in r 2 , where k is the number of €j 
equal to —1. Therefore M is a power of r 2 times a polynomial in r 2 of degree 

fc=0 ^ ' 
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which simplifies to A n . We can factor out the unwanted power of r 2 because 
it was not needed to make equation (|1(J|) hold. 

For semicyclic polygons, the signed sum of the 8j is an odd multiple of 
7r/2. Equation (fTUj) therefore becomes 

II r " cos (^i + e 2# 2 + • • • + t n e n ) = o, 
ej =±l 

which expands to a polynomial relation M'(r 2 , y 2 ,..., y 2 ) = 0. Using series 
expansion again, one finds that M' is monic of degree n2 n ~ 2 in r 2 , and its 
lowest nonzero term has the same degree as that of M. Hence M' is a power 
of r 2 times a polynomial whose degree in r 2 is A n + 2 n ~ 2 = A' n . 

6 Specializations 

Corollaries|Sland|ni which relate the generalized Heron polynomials a n and a' n 
to the polynomials P n (z) and -P^(^), allow us to understand and factor certain 
specializations of a n and a' n . With a little extra work one can describe some 
of the factors explicitly. In this section we offer two such results concerning 
cyclic ra-gons with n odd. 

Let n = 2m + 1 > 5, and consider the constant term of a n \ that is, let 
U2 = 0. Then P n (z) has m — 1 double roots if and only if either (P n \ U2 =o) I ' z 
has m — 1 double roots, or « 3 = and (P n \ U2=U3= o) / z 2 has m — 2 double 
roots. Geometrically, the projective variety 

X — { [v,2 ■ ■ ■ ■ ■ u m : t m+ i : • ■ ■ : t2m+i] | 

P n (z) factors as (b + biz)(c + c\z H h c m _i2 m_1 ) 2 } 

intersects the hyperplane {112 = 0} in two irreducible components, one corre- 
sponding to b = and one corresponding to cq = 0. The second component 
has intersection multiplicity two because X is tangent to {112 = 0} along it. 
Chasing through the geometric interpretation of a n (after Corollary EJ) , we 
find that a n \ U2=0 , as a polynomial in a±, . . . , a n , is an irreducible polyno- 
mial times the square of another irreducible. The factors are not necessarily 
irreducible as polynomials in the side lengths a^, however. 

Proposition 8. If n is odd, the constant term of a n factors as 

OL n 1 16^ 2 =0 = 7n H( a l ± d 2 ± • • • ± O n ) 
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where the product is over all 2 n 1 sign patterns. 

Proof. Heron's formula takes care of the case n = 3, so we may assume n > 5 
and apply the analysis above. By Corollary cyclic n-gons satisfy 

P n {z) = {\-r 2 z){D{r 2 z)/z)\ 

so the factor 7 2 corresponds to = 0, and the other factor corresponds 
to [7 : — r 2 ] = [0:1] and represents projective solutions at r 2 = 00. The 
presence of the linear factors a% ± a-i ± • • • ± a n in the constant term was 
proved in [S], and they correspond to solutions with r 2 = 00: As a signed 
sum of edge lengths approaches zero, the polygon can degenerate to a chain of 
collinear line segments, which has zero area and infinite circumradius. (One 
can easily construct a curve of solutions to the equations (j3J) tending to any 
such point at infinity) For n > 3, the product of these 2 n_1 linear factors is 
symmetric in the a 2 , so by irreducibility, no other factors can appear. □ 

The same kind of analysis applies to a n when the side length a n goes to 
zero, and so t n = a n = 0. It's geometrically clear that the result should be 
divisible by a 2 _x (as the n th side shrinks to zero, it can go either "forward" 
or "backward"), and the algebra confirms it. The intersection of X with the 
hyperplane t n — includes a component of multiplicity two where 

tn-l = ~°~2m + \t^n — 

and P n (z), considered as degree 2m — 3, has m — 2 double roots. Substituting 
the solutions t m = ±2y/a2 m back into the definitions of t m+ \ through ^m-i, 
we recover the definitions of the tj for n = 2m and e = ±1 (see Corollary (SJ) , 
and observe that u m becomes t m . Thus P n {z) specializes to P n _i(z), and so 
a n \a n =o is divisible by (/^-i) 2 ^*^) 2 = a\_ x . 

The other component of X fl {t n = 0} corresponds to solutions in which 
the the leading coefficient r 2 of the linear factor (| — r 2 z) vanishes. We can 
describe these solutions explicitly. 

Proposition 9. If n > 5 and n is odd, then 

«nL=o = «Li U( 16R2 + (o? ± a? ± • • • ± a 2 ^) 2 ), 
the product taken over all sign patterns with (n — l)/2 minus signs. 
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Proof. It suffices to show that all the factors in the product are present; the 
result then follows by comparing degrees (both sides being monic in 16K 2 ). 

Fix generic positive real values for ai, . . . , a n _i and signs e\, . . . , e n _i G 
{ — 1, +1} with €\ = +1 and e j = 0- Recall that cyclic polygons satisfy 

q] + (a 2 /r 2 - 2) Qj + 1=0, 1 < j < n, 

and q±q 2 • • • q n — lj and any solution to these n+1 equations also satifies a n 
with the squared area given by equation (j3J). For sufficiently small positive 
values of a n , the method of §5 shows that there exist solutions (r 2 , qi, . . . , q n ) 
with r 2 < and qj ~ (a 2 /r 2 ) ej for 1 < j < n; furthermore a 2 /r 2 tends to a 
constant. Equation (HJ) now implies 

lim 16K 2 = - 

so the set of solutions with a 2 = r 2 = includes all points on the hypersurface 

i6^ 2 +(E;ri^ 2 ) 2 = - □ 

7 Conclusions 

The generalized Heron polynomials a n for cyclic n-gons and a' n for semicyclic 
(n + l)-gons are defined implicitly by n + 2 equations in n + 1 unknowns: 
n equations © relating the side lengths to the vertex quotients q±, . . . , q n 
and the radius r; equation (jlj) which expresses the squared area K 2 , or equiv- 
alently u 2 = —4K 2 , in the same way; and the equation q x q 2 ■ • -q n — 6. Our 
analysis eliminates the variables qj (and in the cyclic case, r also) at the cost 
of introducing [(n — 5)/2j unwanted quantities u 3 , . . . , u m . The reduction 
in the number of auxiliary variables allows us, for small n, to eliminate them 
by ad hoc means and obtain formulas for ot n and a' n . 

The quantities w 2 , . . . , u m appear on equal footing in this analysis, so 
we could equally well eliminate all but for some k > 2 and obtain a 
polynomial relation, presumably of total degree kA n or kA' n , between Uk 
and the squares of the side lengths. Unfortunately, we do not yet have a 
geometric interpretation for u 3 or the higher u^. 

For large n the goal of eliminating w 3 , . . . , u m seems rather distant, but 
Corollaries El and El still illuminate aspects of the polynomials a n and a' n . 
In particular, Corollary establishes a close relationship between a 2m +i and 
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a 2m+2 , generalizing those between Heron's and Brahmagupta's formulas and 
between Robbins' pentagon and hexagon formulas. 

It may be of some interest to know how our main results were obtained. 
Robbins solved many combinatorial and algebraic problems in his lifetime by 
what he called the "Euler method" : calculate examples, using a computer if 
convenient; discover a general pattern; and prove it, with hints from further 
calculations if necessary. His work on generalized Heron polynomials took 
this approach but was somewhat frustrated by lack of data from which to 
generalize. As explained in |7j, Robbins first found and the closely related 
polynomial a§ by interpolating from several dozen numerical examples, and 
he rewrote them concisely in terms of variables t{ (slightly different from ours) 
by interacting with a computer algebra system. Neither step is particularly 
feasible for a 7 , whose expansion in terms of the symmetric functions has 
almost a million coefficients. It was possible, however, to evaluate certain 
specializations of a 7 and as by interpolation, and to conjecture Propositions 
IH1 and EH For instance, we discovered that the constant term of a 7 (the 
specialization w 2 = 0) is divisible by the square of the discriminant of 

£4 + hz + t 6 z 2 + t 7 z 3 , 

where the ti are defined as in Corollary |S] but with w 2 = W3 = 0. Unfortu- 
nately, the hidden presence of u 3 made it difficult to guess the rest of a 7 . 

We introduced semicyclic polygons and their area polynomials a' n in an 
effort to obtain more data to study. In particular, we noticed that the myste- 
rious cubic discriminant so prominent in appeared already in the simpler 
polynomial a' 3 , and we hoped that whatever new phenomenon arose in a 7 
would also appear in a' 5 , which we could compute by interpolation. In fact 
a 7 turns out to be rather different from a' 5 , but it was the struggle to simplify 
a' 5 that led us to manipulate the relations between the <jj and Tj by hand 
and thence to discover the main identity (Theorem HJ). In the end, the most 
crucial calculations turned out to be those we did on the blackboard. 

Postscript 

David Robbins (1942-2003) had an exceptional ability to see and communi- 
cate the simple essence of complicated mathematical issues, and to discover 
elegant new results about seemingly well-understood problems. He taught 
and inspired a long sequence of younger mathematicians including the two 
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surviving authors. His interest in cyclic polygons began at age 13 when he 
derived a version of Heron's formula. In the early 1990's he discovered the 
area formulas for cyclic pentagons and hexagons. When diagnosed with a 
terminal illness in the spring of 2003, he chose to work on this topic once 
again. Sadly, he did not live to see the discovery of the main identity or the 
heptagon formula. This paper is dedicated to the memory of our friend and 
colleague, whose loss is keenly felt. 
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